Sand stirred by chaotic advection 
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We study the spatial structure of a granular material, A'' particles subject to inelastic mutual 
collisions, when it is stirred by a bidimensional smooth chaotic flow. A simple dynamical model is 
introduced where four different time scales are explicitly considered: i) the Stokes time, accounting 
for the inertia of the particles, ii) the mean collision time among the grains, iii) the typical time 
scale of the flow, and iv) the inverse of the Lyapunov exponent of the chaotic flow, which gives a 
typical time for the separation of two initially close parcels of fluid. Depending on the relative values 
of these different times a complex scenario appears for the long-time steady spatial distribution of 
particles, where clusters of particles may or not appear. 

PACS numbers: 47.52.+j,45.70.Mg 



Many different physical processes can be studied in 
the general framework of transport of finite-size particles 
by an external flow. Sand suspended on the surface of a 
river or driven by the wind, chemical pollutants advected 
by atmospheric flows, and impurities driven by oceanic 
currents, are just but a few of the many important real 
examples that fit into this. 

Recently, much effort coming from the Chaotic Advec- 
tion 1^ and Turbulence ^ communities has been devoted 
to the study of advection of non-reacting inertial particles 
by chaotic or turbulent flows. Though some works con- 
sidered the effect of elastic collisions, see for example , 
the obvious fact that in real systems collisions among 
particles dissipate kinetic energy into heat {inelastic col- 
lisions) has not been, to our knowledge, much subject of 
study. From another perspective, the so-called Granular 
Media community has also widely considered the prop- 
erties of particles colliding inelastically, and externally 
driven, for example, by a shear flow or by vibrating pe- 
riodically the container of the particles . However, the 
influence of an external chaotic or turbulent flow seems 
to be overcome in the literature (see, however, [|j for the 
case of densely packed granular systems). 

In this paper, we try to investigate numerically the 
full problem. Equivalently, we study the effect of col- 
lisions on inertial particles immersed in a flow showing 
chaotic advection 0, or rather, the properties of a di- 
lute granular system (in the literature the term granular 
gas 1^ is widely used) under chaotic stirring. It is by now 
clear that transport behavior in the so-called Batcherlor's 
regime, that is, in the range of scales between the smallest 
typical length scale of the velocity field and the charac- 
teristic diffusion length scale, is equivalent in a turbulent 
flow and a chaotic advection flow. 

We introduce a simple model of granular material con- 
sisting of many particles subject to inelastic mutual colli- 
sions, and immersed in a smooth two-dimensional chaotic 
flow without gravity. Granular particles are assumed to 
be much heavier than fluid. In particular, we focus on 



the steady spatial structure of the system. It is known, in 
fact, that the presence of inertia can induce preferential 
concentration [Q, ^. One may ask what happens if the 
diffusive role of collisions (in the elastic limit) among par- 
ticles is taken into account. At the same time is also well 
known that dissipative collisions originate particle clus- 
tering Therefore, in spite of its simplicity, we will see 
that the model shows many transitions from clustering 
to unclustering depending on the interplay between the 
inertia of the particles, the chaotic flow and the collisions. 

We consider N identical particles of mass to = 1 on a 
two-dimensional domain L x L (with periodic boundary 
conditions) driven by an external velocity field u(x, t) 



dvi(t) 1 , , . , , , 
^^ = --(v,(t)-u(x,(t),i)) 



dx^{t) 
dt 



(la) 
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where i = 1, N, is the velocity of the particle i, r is 
the standard Stokes time, that depends on the particles 
diameter, viscosity, and densities of fluid and particles. 
Hence the term in the r.h.s of Eq. (|la| ) is simply a vis- 
cous Stokes drag. In addition the particles are subject to 
binary instantaneous inelastic collisions according to the 
rule 
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-((v»0-) -Vj(,))-n)n 
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where i and j are the indexes of colliding particles, the 
primes denote the velocities after collision and r G [0, 1] 
is the restitution coefficient (r = 1 for the elastic case). 
This rule ensures the conservation of momentum, while 
the kinetic energy is dissipated if r < 1. 

Some further clarifications on the model are worth 
mentioning. In the absence of collisions the equations in 
(|l|) are simply the equations of motion of a rigid sphere 
in a flow where the Faxen corrections, the added mass 
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term, and the Bernoulli term are neglected This is 
a consistent approximation when the particles of granular 
material are much heavier than that of the fluid, which 
is the case we are considering. Moreover, we adopt the 
Direct Simulation Monte Carlo (DSMC) scheme to inte- 
grate the dynamics of the system The DSMC al- 
gorithm, widely used in numerical hydrodynamics, con- 
sists of a time discretization with fixed time step: in 
every time step there is a free flow-step and a collision 
step. In the free-flow step the motion of particles is in- 
tegrated according to equations without taking into 
account the possible collisions. In the collision step ev- 
ery particle has a probability of colliding with neighbor- 
ing particles (in a disk with a radius smaller than the 
mean free path): the probability of particle i to collide 
with particle j is proportional to the relative velocity 
|vi — Vj|. This scheme has been proved to converge to the 
solution of the Boltzmann equation of the correspond- 
ing hard disk gas and its key feature is (through the 
randomization of collisions) the assumption of Molecu- 
lar Chaos, i.e. lack of correlations for colliding particles, 
P{vi,Vj,t) = P{vi,t)P{vj,t). It is important to note 
that this assumption does not rule out correlations at 
scales larger than the scale of a diameter, and the DSMC 
is often used to study vortices [|2| and other instabili- 
ties. In the context of inelastic gases it has been pre- 
viously used also for the investigation of the clustering 
phenomenon Q. 

A statistically stationary state is reached when the dis- 
sipation of energy because of the inelastic collisions is bal- 
anced with the continuous energy injection coming from 
the external flow |l3) . We determine that the statistically 
stationary state is reached when the typical fluctuations 
of the total energy of the system, E = {vf) (the so-called 
granular temperature), are rather small, typically some 
small percentage of E. As already mentioned, the prop- 
erties of this state are the central objective of this work, 
with special emphasis to the clusterization of particles. 

The formal introduction of the collision operator and 
adimensionalization of the dynamical system Eqs. ( p^ , p^ ) 
will help in order to see the relevance of the distinct 
terms. If we denote T as the typical time of the flow, 
L the typical length and Fcoi/''c the collision operator, 
with Tc the mean collision time, we can rewrite Eq. ( p^ ) 
as follows 

= - v.-u +-Fcoi, 3 
T dt Tc 

where i = t/T, = v,T/L, and u = uT/L. 

For our numerical studies we use the time periodic flow 
u(x, t) = {ua:{x, y, t), Uy{x, y, t)) given by: 

= 2tt La /Tsi\i{2TTy/L), 
Uy = 0, (4) 



if t mod T < T/2, and 
Ux = 0, 

Uy = -2TTLa/Ts\i\{2'Kx/L), (5) 

if t mod T > T/2. 

Periodic boundary conditions are assumed, and for a > 
0.4 a unique chaotic region without KAM tori (numer- 
ically distinguishable) is obtained iQ. The Lyapunov 
exponent for the flow is given by A = 1.96 ln(3.35Q;)/r 
and introduces a new time scale, Tf — 1/A, for the ex- 
ponential separation rate of close fluid parcels. Our re- 
sults has been checked to be independent of the kind of 
flow. In particular, we have alternatively used the cel- 
lular flow derived from the stream function ^'(a;,y,t) — 
U g^j^^i "^^+-^'^°'^('^*) ^ gjjj^npt)^ where U is the maximal ve- 
locity of the flow, n the number of cells, B the strength 
of the time dependent forcing and w its frequency (l^. 

Depending on the relative values of the above intro- 
duced time scales: r, T, r/ and Tc, a very complex sce- 
nario appears for the spatial distribution of the particles 
in the stationary state. Some preliminary general results 
can be anticipated. For example, it seems obvious that 
when Tf < Tc the exponential separation produced by the 
chaotic driving of the fluid may avoid any clusterization. 
The opposite may happen if Tc < Tf as the mechanism 
of aggregation due to inelastic collisions resists the dis- 
persing effect of the chaotic flow. However, as we will see 
it is also very important to take into account the com- 
paction due to the inertia. In any case, discussing all 
the different possibilities (taking also into account some 
values of the parameter r) goes beyond the scope of this 
Paper. Thus, we just restrict ourselves to study some of 
the most interesting cases. 

First we consider the case of irrelevant inertia, t «T . 
Obviously, in the absence of collisions, particles fol- 
low strictly the fluid and no aggregation appears (see 
Fig. |l|a)). In fact, r is the relaxation time of the particles 
to the driving flow; without collisions the particles are 
continuously flowing through the entire spatial domain 
since no KAM tori are present in the flow. However, 
everything changes when the collisions are taken into ac- 
count. Thus, if Tc < T, the last term in Eq. ^) is the 
most relevant one. If we have also that Tc < Tf an aggre- 
gation of the particles is always observed for any value 
of r (even, and most noticeably, at r = 1). This can be 
understood from Eq. (^), since now the velocity of the 
particles can be written as: 

v«u+-Fcoi. (6) 

Tc 

Thus, the effect of collisions is similar to that of iner- 
tia when considered alone: the velocity of the particles 
is modified respect to that of the fluid parcels, then it is 
no longer compressible (V ■ v ^ 0), and there are re- 
gions where particles aggregate. Figs. |^b) and ^) show 
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FIG. 1: Spatial distribution of particles in the steady state, 
when the inertia is irrelevant (or almost), r < T. a) Case 
with no collisions. T = 10, a — 10 {jf ~ 2), and r = 0.1. 
The number of particles is A*' = 500 and the final time 500. b) 
Same as before but with almost elastic collisions and « 
Tf. Here r = 0.99, Tc = 0.01, r/ ~ 2. c) Same as before but 
with inelastic collisions, r — 0.6. d) Also the case of inelastic 
collisions, but now r/ << Tc- The number of particles is 
iV = 200 and T = 0.2, a = 1 {rf ^ 0.08), r = 0.2, = 0.1, 
and r — 0.1. Final time also 500 units. 



a pattern of the distribution of particles for r = 0.99 
and r = 0.6, respectively. It is also very important to 
note here a specific behavior of elastic or quasi-elastic 
collisions (r « 1). The effect of multiple elastic col- 
lisions is equivalent to macroscopic diffusion, which, if 
large enough, may induce a dispersing mechanism that 
prevents clustering. Summing up, one can say that elas- 
tic collisions induce a clustering inertia- like effect, that 
disappears when diffusion is strong enough. Quantita- 
tively this is controlled by the adimensional parameter 
t/tc since the macroscopic diffusivity D oc 1/tc. Increas- 
ing t/tc clustering like that of Fig. |l]b) disappears and 
a homogeneous distribution is attained (see discussion 
below on the P(m) function and Fig. |^)). 

On the contrary, when Tf < Tc (and also Tc << t) the 
strong chaoticity of the flow dominates and avoids the 
formation of any clustering in the system, see Fig. |l|d). 
Finally, it is obvious from Eq.(||) that the limit t << Tc 
is equivalent to the absence of collisions, thus no cluster- 
ization appears. 

A more quantitative analysis of clustering follows. This 
is performed via a particle- in-cell histogram, i.e., after di- 
viding the system in M small boxes (we use N/M = 5) 
an histogram is made with the number of boxes contain- 
ing TO particles, denoted the corresponding function as 
P{m). For an homogeneous system of particles P{m) is 
a Poisson distribution ea;p(— A)A'"/m!, where A = N/M . 
As the clustering is stronger developed, the deviations 




b) 



O no collisions 

□ r=0.99(T^<t) 

o r=0.6 (t_,<t:) 

< r=0.99 (T_,«T) 

A Strong chaos (r=0. 1 ) 



o no collisions 

□ r=0.99 

o r=0.75 

A Strong chaos (r=0.75) 




'OoOOOgOo 



o^oo! 



FIG. 2: Distribution of clusters for the same cases shown in 
figure 1 (above) and 3 (below). Solid line is the analytical 
Poisson distribution. 



from a Poissonian are more evident. Fig. [^a) calculates 
P{m) for the patterns shown in Fig. ^ and one more case. 
We observe the appearance of clustering for r < 1, and 
even for r « 1. Also, it is shown how for the elastic case 
the Poisson distribution is approached when t/tc << 1. 




FIG. 3: Spatial distribution of particles in the steady state, 
when the inertia is relevant, r > T. a) Corresponds to the 
case without collisions. Here T — 1, a — 0.5, and r = 10. 
The number of particles is A*' = 500 and the final time 500. 
b) Same as before but now with almost elastic collisions, with 
Tc « Tf. Here r = 0.99, Tc = 0.1, r/ ~ 1. c) Same as before 
but with inelastic collisions, r = 0.75 (here the final time is 
175 as the program breaks because of collapse earlier, d) Also 
the case of inelastic collisions, but now r/ « Tc. The number 
of particles is = 200 and T = 0.2, q = 1 (r/ f« 0.08), r = 2, 
Tc — 0.2, and r = 0.75. Final time also 500 units. 

It is very relevant to know the effect of collisions in the 
context of transport of inertial particles in chaotic flows. 
Therefore, we next consider the general case where the 



4 



inertial term is relevant, t > T (but otherwise r small 
enough such that the l.h.s is not negligible with respect to 
first term in the r.h.s). It is well-known |^ that, in the ab- 
sence of collisions, heavy particles accumulate in regions 
of low strain and high vorticity, the so-called preferential 
concentration phenomenon. In our case, the clustering 
steady pattern is shown in Fig. ||a) . After including col- 
lisions one has: if r << Tc no effect is observed as the 
last term in the r.h.s of Eq. (||) can be neglected. A 
different behavior is observed when Tc << t. Regardless 
of the value of r/, in the limit of elastic collisions (r « 1) 
no aggregation of the particles is observed. Fig. ||b). As 
explained before, when t/tc is large, the continuous elas- 
tic colliding processes of the particles induce a diffusive 
effect that prevents clustering. On the contrary, as we go 
into the inelastic regime, r < 1, the chaotic time scale, 
T/, now plays a very important role. Thus, when < Tf 
particles aggregate although in a very different pattern 
to that obtained in the absence of collisions (Fig. ^)). 
In particular, at variance with clustering in the absence 
of colhsions (see Fig. ^)), clusters are created and de- 
stroyed continuously, and moving with the flow. In this 
case, the particles aggregate thanks to the inelastic colli- 
sions and because the flow is not chaotic enough to dis- 
perse completely the clusters. The opposite is observed. 
Fig. ||d), when Tf < Tc, the clusters do not resist the 
dispersion of the flow and the density of particles is ho- 
mogeneous. The corresponding P(m) function is plotted 
in Fig. |b). 

For completeness, we briefly mention the limit t » T . 
This is completely equivalent to a freely cooling granular 
gas since the external flow is irrelevant. Thus, clustering 
is usually observed for inelastic collisions p^ . 

In this letter we have studied numerically the station- 
ary spatial structure of a granular material advected by a 
chaotic flow. We have seen that the relative importance 
of collisions, inertia and chaoticity of the flow gives rise 
to transitions from homogeneous to inhomogeneous den- 
sity of grains. We think that our simple model is relevant 
for real physical phenomena that involve the transport of 
a large amount of heavy particles, like for example those 
appearing in many geophysical situations. 
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